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Abstract 



The synchrotron sclf-Compton (SSC) models and External Compton (EC) models of AGN jets with 
continually longitudinal and transverse bulk velocity structures are constructed. The observed spectra 
show complex and interesting patterns in different velocity structures and viewing angles. These models 
\ are used to calculate the synchrotron and inverse Compton spectra of two typical BL Lac objects (BLO) 

(Mrk 421 and 0716+714) and one Flat Spectrum Radio Quasars (FSRQs) (3c 279), and to discuss the 
implications of jet bulk velocity structures in unification of the BLO and FR I radio galaxies (FRI). By 
calculating the synchrotron spectra and SSC spectra of BL Lac object jets with continually bulk velocity 
structures, we find that the spectra are much different from ones in jets with uniform velocity structure 
under the increase of viewing angles. The unification of BLO and FRI is less constrained by viewing angles 
and would be imprinted by velocity structures intrinsic to the jet themselves. By considering the jets 
with bulk velocity structures constrained by apparent speed, we discuss the velocity structures imprinted 
on the observed spectra for different viewing angles. We find that the spectra are greatly impacted by 
longitudinal velocity structures, becasuc the volume elements arc compressed or expanded. Finally, we 
present the EC spectra of FSRQs and FR II radio galaxies (FRII) and find that they are weakly affected 
by velocity structures compared to synchrotron and SSC spectra. 
, Key words: galaxies: jets - radiation mechanisms: nonthermal 
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' 1. Introduction bulk velocity along the jet axis. Their model is also 

0^ . used to explain the X-ray phenomena in kiloparsec jets 

It is well known that the AGN jets radiate the strong of powerful radio galaxies and quasars (Georganopoulos 
nonthermal emission from the radio to the gamma ray & Kazanas 2003b). Girolctti et al. (2004a) have pro- 
v~j ■ range. Their spectral energy distribution (SED) con- posed a transverse velocity structure model, described 
rS ' sists of two bumps, attributed to the synchrotron and by a slower external flow surrounding a faster spine, to 
^ ■ the inverse Compton (IC) emission of ultrarelativistic explain a limb brightening morphology in Mrk 501 ob- 
■ ' particles. The simple one-zone homogeneous jet model served by VLBI. Ghiscllini & Tavccchio (2008) have used 
usually provides a good framework to explain the emis- a needle/jet model, which harbors small active regions 
sive spectra of AGN jets. However, the realistic case (needls) inside a large jet, to account for the 2-3 min- 
is that the jet has the bulk velocity structures in longi- utes fast variability of PKS 2155-304. Recent observation 
tudinal and transverse directions to its axis (Blandford of TcV blazars raise new questions regarding the one-zone 
& Lcvinson 1995). To explain the spectral properties uniform jet model. Such as, the unrcally high Lorentz fac- 
of some AGNs several authors have proposed many in- tor (about 50) demanded by avoiding 7-7 absorption for 
homogeneous jet models (Marscher 1980; Ghiscllini et TeV blazars (Krawczynski et al. 2002; Konopelko et al. 
al. 1985; Georganopoulos & Marscher 1998; Wang et al. 2003; Henri & Sauge 2006; Begelman et al. 2008; Finke 
2004; Georganopoulos et al. 2005; Tavecchio & Ghiscllini et al. 2008) can not match the much slower speeds re- 
2008). For example, to reproduce the spectra of TeV quired by Very Long Baseline Interferometry (VLBI) ob- 
blazars, Georganopoulos & Kazanas (2003a) have pre- servations (Pincr & Edwards 2004; Piner ct al. 2008). 
sented a jet model characterized by the deceleration of Therefore, it is argued that the jet should have veloc- 
ity structures to mediate the inconsistence of bulk mo- 
tions obtained by different 
2000; Trussoni et al. 2003). 
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For unification of BL Lac objects and FR I radio galax- 
ies, the velocity structures have also been suggested by 
Chiabcrge et al. (2000). By analyzing the core emis- 
sions of FR I radio galaxies observed by Hubble Space 
Telescope, they have found that the observed fluxes of 
FR I nucleus are over- luminous by a factor of 10 — 10 4 
than ones predicted by simple one-zone model, in radio 
and optical bands. Therefore they proposed a jet velocity 
structure, which a fast spine is surrounded by a slow layer, 
to reconcile the unification scheme. 

In this paper we present the SSC and EC models of 
AGN jets with continually longitudinal and transverse 
bulk velocity structures. Our study is not completely new, 
but we include new ingredients and results. Firstly, the 
observed apparent speed is obtained by flux convolution 
over the whole emitting region with bulk velocity struc- 
tures. Secondly, each part of the jet contributes its spec- 
trum to the observed SED, some important characteris- 
tics are ignored in simple models with two-zone velocity 
structures compared to the model with continual veloc- 
ity structures. Thirdly, some parameters of jet velocity 
structures can be constrained by the jet power or appar- 
ent speed , our models are more realistic to discuss the 
emissive properties of the jet with velocity structures than 
previous models. Fourthly, the jet velocity structures can 
reasonably solve the problem of the unification of BL Lacs 
and FR I faced by all one-zone models. Fifthly, we find 
that the effects of velocity structures upon the observed 
synchrotron and SSC spectra strongly depend on the view- 
ing angles. Sixthly, we find that the longitudinal velocity 
structures, as the volume elements are compressed or ex- 
panded, have great influences upon the observed spectra. 
Finally, the observed EC spectra are weakly affected by 
velocity structures compared to the observed synchrotron 
or SSC ones. 

In Section 2 we present the velocity structures of the 
jet for the models calculation and we place these modifi- 
cations in the context of emission mechanisms in Section 
3. We show the model application to two BL Lac objetcs 
(Mrk 421 and 0716+714) controlled by the same jet power 
in Section 4, focusing on the unification of BLO with FRI 
sources. We then discuss the apparent speeds that govern 
the velocity structures of the jet to constrain the observed 
spectra in Section 5. In Section 6, we employ the SSC and 
EC models to one quasar (3c 279), and discuss the effects 
of the velocity structures upon the EC spectra. We finish 
with discussions and conclusions in Section 7. 

2. Velocity Structures 

We only select a part of the stationary jet to study. For 
simplicity, we adopt the following assumptions: 1. The 
length of the studied region is a small scale (i.e., about 
1 x 10 17 cm) compared with the whole jet. Although the 
whole jet may be conical or other morphologies, a cylinder 
jet is a good approximation within a short length scale. 
2. In the studied region, the jet is assumed to have the 
bulk velocity structures in the longitudinal or transverse 
directions to the jet axis. The velocity structures are 



not resolved by the VLBI observation. 3. We assume 
that the jet is stationary and the relativistic electrons 
are injected continuously with a power-law energy spectra 
A( 7 ,z ) =N Ql - s exp{-^ l; ) (Finke et al. 2008). Here we 
mainly pay attention to the influence of the velocity struc- 
tures on radiative spectra, we ignore the radiative energy 
loss of the electrons in the studied jet scale. However, 
for the jet with the longitudinal velocity structures, as its 
volume elements are compressed or expanded, the elec- 
tron energy distribution and magnetic field in the jet will 
be changed follow as Georganopoulos & Kazanas (2004): 
7 (z)= 70 (z/z )-?/ 3 , JV( 7 ,z) = JV( 7 ,z )(z/zor« (s+2)/3 and 
B(z) = Bo(z/zo) where the £ is the index of longitudi- 
nal velocity structures (see the next part), the subscript 
"0" denote the values at the studied jet base. 4. For trans- 
verse velocity structures, the magnetic fields are assumed 
to be a isotropic distribution in the comoving frame. 

2.1. Longitudinal Velocity Structures 

For a jet with longitudinal velocity structures, we take 
the Lorentz factor along the jet to be T(z) = To(-^)*, 
where To is the Lorentz factor at the base, the £ is the 
index of longitudinal velocity structures. This power-law 
form was earlier given by Marscher (1980), subsequently 
adopted by Ghisellini et al. (1985); Georganopoulos & 
Marscher (1998) and Li & Wang (2004). £ = 0, £ > 0, 
and £ < correspond to transition, acceleration, and de- 
celeration phases respectively. 

VLBI measurements of apparent motion for the parsec- 
scale radio knots have often been employed to con- 
strain a combination of Lorentz factor and viewing an- 
gles (Vermeulen & Cohen 1994). Although VLBI moni- 
toring of the radio knots in blazar jets has revealed several 
sources containing knots with apparent speeds extending 
out to 30c (Piner et al. 2007), the typical values for TeV 
blazars are found to be much more modest, under 5c (e.g. 
Giroletti et al. (2004b); Piner & Edwards (2004)). Gopal- 
Krishna et al. (2004) have showed that the slow apparent 
speeds of the knots of blazars observed by VLBI can be 
reconciled with the extremely relativistic bulk motion in- 
ferred from TeV flux variations (Krawczynski et al. 2002), 
if one considers a modest full opening angle for the parsec- 
scale jets. They also have showed that the actual viewing 
angles, 8, of such conical jets from the line-of-sight can 
be substantially larger than those commonly inferred by 
VLBI proper motion data (Gopal-Krishna et al. 2006). 
Recently, they have evaluated the role of the jet open- 
ing angle on certain key parameters (i.e. the viewing an- 
gle, the apparent speed and Doppler factor) that are in- 
ferred from VLBI radio observations of blazar nuclear jets 
(Gopal-Krishna et al. 2007). In these papers, they have 
argued that the Doppler boosting of an ultrarelativistic 
jet, as well as the apparent proper motions, can greatly 
vary across the jets cross-section, and then it is important 
to carry out an integration of various quantities across the 
jet cross-section. Therefore, they have performed an inte- 
gration of the (boosted) flux-weighted apparent velocity 
over the jet cross section to obtain the weighted observed 
value of the apparent velocity of the jet, when the width 
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of a knot cannot be resolved by VLBI observations. 

For longitudinal velocity structures the observed appar- 
ent speed (Papp = i-pcose ) ^ s decided by (Gopal-Krishna 
et al. 2004; Gopal-Krishna et al. 2006; Gopal-Krishna et 
al. 2007) 
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where the F sync is the observed flux, and the dF sync (z) 
is the comoving flux in z. In this equation, we take the 
radio spectral index to be 0. 

Through the radio observation, Giroletti et al. (2004b) 
have found that the mean value of T is 3 in the TeV 
sources. Piner & Edwards (2004) have presented that 
the jets of TeV blazars have only mildly relativistic mo- 
tion which is less than 5c on parsec-scalcs. Pincr et al. 
(2007) have found that the apparent speed distribution 
shows a peak at low speed for BL Lac objects. We then 
take the speed of (3 a pp.obs = 3.6c to limit Tq. 

2.2. Transverse Velocity Structures 

For a jet with transverse velocity structures, we suppose 
the jet bulk velocity to be T{R) =T C - (T c - T m )(R/R )^, 
where F c and r m are the Lorentz factor of the central and 
marginal region respectively. In this work, we only adopt 
£ to be and 1 , corresponding to different transverse ve- 
locity structures. 

We take the apparent speed given by (3 a pp,obs = 3.6c 
to restrict T c and T m . We also use the jet power given 
by (Celotti et al. 1997; Ghisellini et al. 2005) P jet = 

J^°r 2 (R)UP(R)c2wRdR to limit T c and T m , where U = 
Ub + U e + Up is the total energy density in the jet frame. 

3. Radiation Mechanisms 

In this section, the radiation models of continually 
longitudinal and transverse velocity structures are con- 
structed respectively 

3.1. Radiation Model of Longitudinal Velocity Structures 

We present a modified homogeneous SSC model to re- 
produce the entire spectral energy distribution. In the 
following, a prime accent denotes a parameter expressed 
in the comoving frame. 

For an isotropic electron distribution, the synchrotron 
emission coefficient is given by 



(2) 



where v s is the frequency in the comoving frame, and 
P e (v s , r f) is the mean emission coefficient for a single elec- 
tron averaged over an isotropic distribution of pitch an- 
gles. The differential synchrotron luminosity between z 
and z + dz is then given by 



dL s (v s ,z,8) = 4ir 2 R 2 j' s (v s ,z)e 



'dz, 



(3) 



where Rq is the radius of jet, a s {v s ) is the absorption 
coefficient (Rybicki & Lightman 1979) as follows 



2 d , N e(l,z) , 

8nm e v' 2 J 7 ^^^ P ^,7)d7.(4) 

The Doppler factor is given by 6(z, 9) = r {z) [ i-l(z)cos6 ] an d 
r(z) = r (^)«, where T(z) = [1 - ^(z)]" 1 / 2 is the bulk 
Lorentz factor and (3(z) is the velocity in the unit of c. 
Then the synchrotron flux density in the observe frame is 
simply given by 



dF. 



sync 



( Vs ,z 1 9) = -\ 1 S^ +a \z,6)[l + z rs }dL s (iy' s ,z,6),(5) 



And 



where a is the spectral index and is taken as 1.0. di is the 
luminosity distance and z rs is the redshift. We transfer 
v s into the observer's frame following as: 



5(z,9) 
v s [z,0) = -v< 



1 + Z r ,\ 



(0) 



Finally, we get the total flux density through the integra- 
tion of v s {z,9) by z. 

Calculating the differential SSC emission between z 
and z + dz, we need to compute soft photon density pro- 
duced at layers with different Lorentz factors, which is 
given by u ph ~ u ph T 2 ela (Dermer 1995; Georganopoulos 
& Kazanas 2003a). T re i a is the relative Lorentz factor be- 
tween soft photons and IC electrons. u* h is soft photon 
energy density in the local frame. The photon energy den- 
sity u (v g ,z,9) is given by a sum of a local contribution 
(Kataoka et al. 1999): 

Uioc[v S i z ) = i — r-r~' — d 1 ~ ex P\- a s{ v si z )Ro)l ( 7 ) 
4 ca s \y s ,z) 

and an external contribution: 



u ext{ v si z ^) 



j s (^,z)7ri? 2 



2 f-z-R 



^rela( Z S} Z i 



z+Ro 



exp[-a s (u s ,z)(z - z s )]dz s 

(Z-Z s f 

exp[-a s {v' s ,z){z s - z)]dz 
' (z s -z) 2 



+ 



(8) 



Then the emission coefficient is given by (Inoue & 
Takahara 1996; Kataoka et al. 1999; Katarzynski et al. 
2001) 



j c (v c ,z,6) = ^-e c q(e c ,z, 



(9) 



where q(e c ,z,9) is the production rate of the differential 
photon: 

q(e' c ,z,6)=f de' s n ph (e' s ,z,9) J d^N(-f,z)C(€ c7 j,e' s ).(10) 

C(e c ,7,e s ) is the Compton kernel given by Jones (1968) 

r i 27rr 2 c 
C(e c , 7 ,e s ) = _-£_[2/rfn(«) + (1 + 2«)(1 - k)+ 



(4e a 7K) 
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2(i+4 £ ; 7 «) 
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and r e is the classical electron radius. The integration of 
(10) follows the condition as: 



e s < £ c < 7 



l + 4e; 7 ' 



(13) 



Then, we can integrate dL c {v cl z,9) to obtain the total 
flux density. 

In calculating the differential EC emission between z 
and z + dz, we need to know soft photon energy den- 
sity produced by external radiation field, we mainly con- 
sider the photons reprocessed by Broad Line Region 
(BLR)(Sikora et al. 1994). The photons are assumed 
to be distributed as a blackbody with temperature T ext . 
In the comoving frame of each layer the mean soft pho- 
ton energy is blueshiftcd to e ext ~ T(z)e ex t ~ r(z)/csT ea;t , 
where the fee is the Boltzmann constant. The photon 
energy density is given by u ph ext ~ u ph . ext T 2 (z). 

3.2. Radiation Model of Transverse Velocity Structures 

To the transverse velocity structures, we replace the 
equation (3) by 



dL s {v s ,R,6) = 8TT 2 Rj s (v s ) (z max -zo)e 
Soft photon energy density is given by: 

27rj>;) 



s (" s )Bo 



-dR.(14) 



u (v s ,R, 



dz 



dR s T 2 ela (R s ,R,9) 



R s exp[ 



as (Vs)(R 2 + Z 2 ) 1/2 \ 



15) 



/o (R 2 + z 2 ) 

where R s is the radius of producing the soft photons. 
Since the studied length is assumed to be short com- 
pared to the whole jet, we take the soft photon density 
at 2t "°|~ 2q in different z as an approximation . We find 
that this simplicity is reasonable in the following calcula- 
tion. 

4. Velocity Structure Constrained by Jet Power 

In this section, we apply the transverse velocity struc- 
tures constrained by jet power to the SEDs of Mrk 421 and 
0716+714, and discuss their properties in the unification 
of BL Lac objects and FR I radio galaxies. 

We select the well observed Mrk 421 that is a typically 
high energy peaked BL Lac (HBL). The SED of Mrk 421 
is obtained by fitting the observed core data (for the qui- 
escent state of Mrk 421, the data is from the Macomb et 
al. (1996) and NED) with our model. The SEDs of its par- 
ent population are then given by increasing the viewing 
angles. In this work, we adopt two kinds of transverse ve- 
locity structures £ = (uniform velocity structures) and 
C = 1 (bulk Lorentz factor linearly lessen from the cen- 
ter to the margin). The observed data of Mrk 421 are 
fitted using the parameters listed in Table 1. Its corre- 
sponding SEDs in different viewing angles are calculated 
to show the implications of velocity structures. Konopelko 



et al. (2003) have used the one-zone SSC model to fit the 
lower state SEDs of Mrk 421 and obtained the parame- 
ters of s = 1.75, ^max = 3 x 10 5 , B = 0.10G; Blazejowski 
et al. (2005) have presented the parameters of 5 = 10, 
B = 0.405G, R = 0.7 x 10 1( W The parameters in our 
model are comparable to their ones. We set F c = r m = 10 
for the jet with uniform structure and F m = 3 for the jet 
with velocity structure (Ghisellini et al. (2005); Girolctti 
et al. (2004a); Tavecchio & Ghisellini (2008)). Assuming 
the same jet power for two structures, we can obtain F c 
for the jet with velocity structure. The SEDs of Mrk 421 
are shown in Fig. 1. The SEDs of two velocity structures 
show very small difference in the small viewing angle (3°). 
Under the viewing angle of parent population (FR I) (60° 
), the SEDs have large difference, in which the flux with 
velocity structure is much larger than ones with uniform 
structure. The velocity structure can resolve the prob- 
lem that the observed FR I nuclei are over-luminous by a 
factor of 10 — 10 4 than ones predicted by simple one-zone 
model in optical and radio bands. In order to compare 
the observed and predicted fluxes of parent population, 
we also present two vertical lines denoting the optical (V 
band) and 5 GHz radio band respectively in Fig.l. It is 
noted that the radio fluxes given by the model arc extrap- 
olated from the infrared fluxes. 

0716+714 is a typical low energy peaked BL Lac (LBL). 
We also list the parameters in the Table 1 . Giommi et al. 
(1999) gave the magnetic field to be B > 0.9G. Tagliaferri 
et al. (2003) got the parameters of B = 2.5G, R = 2 x 10 16 , 
8 = 3.4°. These parameters are consistent with ones given 
by our model. The SEDs of 0716+714 given by the model 
are shown in Fig. 2 . It is shown that the SEDs have 
obvious difference for two velocity structures under large 
viewing angles (40° or 60°). 

In the Fig. 3, we give the luminosity ratio between the 
jets with C = 1 an d C = in three bands: radio (9GHz), op- 
tical (V band), and X-rays (IkeV) for Mrk 421, 0716+714, 
and their parent population. It is shown that the ratio 
can reach 2 — 9 in three bands under large viewing an- 
gles corresponding to parent population. Understanding 
these properties, we present the variations of Doppler fac- 
tors along the jet radius R for different velocity structures 
under the viewing angles of 3° and 60°, shown in Fig. 4 
and 5. Under 3° (see the Fig. 4), the Doppler factor of 
the jet with velocity structure rapidly decreases with R 
and is smaller than ones of the jet with uniform velocity 
structure when R > 0.6Rq. However, the increase of the 
Doppler factor with R appears under 60° (Fig. 5). 

In the unification scheme, it is believed that the appear- 
ance of AGNs strongly depends on the viewing angles and 
obscuration instead of intrinsic physical properties. For 
low luminosity radio-loud objects, such as BL Lac objects 
and FR I radio galaxies, the relativistic beaming effects 
caused by viewing angles should play an important role in 
observation. However, Chiaberge et al. (2000) have found 
that the observed FR I nuclei are over-luminous by a fac- 
tor of 10 - 10 4 than ones predicted by simple one-zone 
model in the radio and optical bands. They have argued 
that a radial velocity structure, which a fast spine sur- 
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Fig. 1. The approximately fitting for data of Mrk 421 (from 
Macomb et al. (1996) and NED) and the extrapolated SEDs 
of parent population. The solid lines represent the emission of 
jet with uniform velocity structures (£ = 0) and the dash lines 
denote the jet with velocity structures (£ = 1). The upper 
two curves report the emission at angles of 3°, the mid and 
lower two curves are the cases of parent population of BL Lacs 
extrapolating the viewing angles to 40° and 60°. 
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Fig. 3. The ratio of luminosity between f = 1 and £ = under 
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Fig. 2. The approximately fitting for NED data of 0716+714 
and the extrapolated SEDs of parent population. The solid 
lines represent the jet with uniform velocity structures (£ = 0) 
and the dash lines denote the jet with velocity structures 
(£ = 1). The upper two curves report the emission at angles 
of 2°, the mid and lower two curves are the cases of parent 
population of BL Lacs extrapolating the viewing angles to 
40° and 60°. 
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Fig. 4. Under Q = 3°, the Doppler factors vary with R for 
different velocity structures. Solid line denotes the case of 
uniform velocity structures and the dash line is in velocity 
structures. 

rounded by a slow layer, can reconcile the above contra- 
diction. We apply the velocity structures described pre- 
viously to the unification of BLO and FRI. In Fig. 6, we 
show the debeaming trails with different velocity struc- 
tures in the radio-optical luminosity plane for Mrk 421 
and 0716+714. For Mrk 421, when the viewing angle is 
increased to 60°, the predicted optical and radio luminos- 
ity nearly come into the region of FRI in case of velocity 
structure, however they are less-luminous in case of uni- 
form structure. For 0716+714, the predicted luminosity 
fully fall into the region of FRI in different velocity struc- 
tures. 

The a ro -a ox planes for Mrk 421, 0716+714 and FRI 
are presented in the Fig. 7, in which the lines denote the 
debeaming trails of Mrk 421 and 0716+714. The plotted 
data of FRI and BLO are from the Fossati et al. (1998) 
and Trussoni et al. (2003). For 0716+714, the debeamed 
indices fall into the region of FRI under different veloc- 
ity structures. However, the debeamed indices of Mrk421 
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Fig. 5. Under (9 = 60°, the Doppler factors vary with R for 
different velocity structures. Solid line denotes the case of 
uniform velocity structures and the dash line shows the case 
of velocity structures. 



Fig. 7. Debeaming trails in the a ox - a ro plane for Mrk 421 
and 0716+714. The solid denote £ = and dash lines is f = 1. 
The lower pentagrams represent BL Lac, while the upper ones 
describe the parent population of BL Lac observed at 9 = 40° 
and 60°. The data of FR I, HBL, and LBL are given by 
Fossati et al. (1998); Trussoni et al. (2003). 
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Fig. 6. Debeaming trails with the models of continual veloc- 
ity structures in the radio-optical luminosity plane for Mrk 
421 and 0716+714 (the solid lines denote C = and dash ones 
are £ = 1), where upper pentagrams correspond to the lumi- 
nosity of BL Lacs at 9 = 2° or 3° and lower ones are the pre- 
dicted luminosity of parent population at 9 = 40° and 60° . 
The three regions enclosed by curves respectively describe 
the range for three samples, i.e., HBL (dash curves), LBL 
(dot curves), and FR I (solid curves) in realistic observations 
(Chiaberge et al. 2000). 

are marginally in the region of FRI. This implies some 
intrinsic physical difference exist between HBL and LBL. 

5. Velocity Structure Constrained by Apparent 
Speed 

In this section, we use the observed apparent speed to 
constrain the jet velocity structures and discuss its effect 
on the SEDs of BL Lac objects shown in Mrk 421 and 
0716+714. 



We take (3 ap p = 3.6c to limit longitudinal and trans- 
verse velocity structures and calculate the SEDs of the 
jet with velocity structures to fit the observed data of 
Mrk 421 and 0716+714. Under /3 app = 3.6c, adopting 3° 
or 2° as the viewing angles we obtain Tq for three kinds 
of longitudinal velocity structures, e.g. uniform (£ = 0), 
decelcrating(£ = —1), and accelerating (£ = 1), which are 
listed in the figures. The SEDs of Mrk 421, 0716+714 and 
their parent populations under different velocity struc- 
tures are shown in Fig. 8, 9, 10, and 11. It is shown that 
the elementary volume expansion (£ > 0) or compression 
(£ < 0) change the electron energy distribution in longi- 
tudinal bulk velocity structures and cause the parameters 
fitting the observed data to be large difference between 
uniform and velocity structure models (see the Table 1). 
In the decelerating velocity structures, we use smaller pa- 
rameters of No, jmax and B to fit the observed SEDs 
compared to the uniform structure. On the contrary, we 
must use larger parameters in the accelerating models to 
fit the SEDs. When the velocity structures exist, the flux 
vF v decreases slowly along the increase of the viewing an- 
gles (such as 40° or 60°). However, the flux vF v decreases 
quickly in the uniform structure. In the Fig. 12 and 13, we 
give the luminosity ratio between different viewing angles 
in the same band. In the uniform jet, the ratio changes 
greatly with the viewing angles. The energy equiparti- 
tion between electron and magnetic field could not be 
true(Hardcastle et al. 1998). We find that for the accel- 
eration velocity structures U e /UB has larger value (about 
300) at the outside of studied jet scale, however, for the 
deceleration velocity structures the U e is larger than Ub at 
the inner of studied jet scale, along the jet, Ub is became 
dominant at the outer (see the Fig. 14). 

To the transverse velocity structures (Fig. 15 and 16), 
we adopt two kinds of velocity structures, e.g. £ = (uni- 
form velocity structure) and £ = 1 (bulk Lorentz factor 
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Fig. 8. Under f3 app = 3.6c, the approximately fitting for data 
of Mrk 421 (from Macomb et al. (1996) and NED) and the 
extrapolated SEDs of parent population (Note that we use 
different parameters for different velocity structures, see the 
Table 1.). The solid lines represent the emission of jet with 
uniform velocity structures (£ = 0) and the dash lines denote 
the jet with velocity structures (£ = —1). The upper curves 
show the emission at angles of 3°, the red and green curves 
present parent population of BL Lacs with the viewing angles 
of 40° and 60° . 



Fig. 9. Under (i apP = 3.6c, the approximately fitting for data 
of Mrk 421 (from Macomb et al. (1996) and NED) (Note that 
we use different parameters for different velocity structures, 
sec the Table 1.) and the extrapolated SEDs of parent pop- 
ulation. The solid lines represent the emission of jet with 
uniform velocity structures (£ = 0) and the dot lines denote 
the jet with velocity structures (5 = 1). The upper curves 
report the emission at angles of 3°, the red and green curves 
present parent population of BL Lacs with the viewing angles 
of 40° and 60°. 



linearly decreases from the center to the margin). We set 
T m = 3 and obtain r c restricted by (3 apP = 3.6, which is 
listed in corresponding figures. The SEDs of Mrk 421 and 
0716+714 in different viewing angles are shown in Fig. 15 
and 16. We show that the lower luminosity is common to 
the jets with uniform velocity structure in large viewing 
angles, and that the way of solving the paradox of BLO 
and FRI unification is to consider the velocity structure, 
where all photons are not produced in same Doppler fac- 
tor. 

6. EC Spectra with Velocity Structure 

Piner et al. (2007) have proposed that the mean fastest 
apparent speed for the quasars is 6.8 ± 1.1c, as the pre- 
vious discussion we will use this observed apparent speed 
to constrain the parameters of velocity structures (i.e., 
r = 19.6 for the £ = -0.5; T = 8.4 for the £ = 0). In the 
Fig. 17, under (3 app = 6.8c, we present the fit for the data 
of FSRQ 3c 279 (Inoue & Takahara 1996). In the model, 
we set T ext = lOeV (near the energy of hydrogen Lyman- 
a photons) and U ex t = 6.6 x 10 -5 erg cm -3 for external 
photons, other parameters are shown in the Table 1. The 
solid lines represent the jet with uniform velocity struc- 
tures (£ = 0) and the dot lines denote the jet with velocity 
structures (£ = —0.5). The upper curves show the emission 
with the angle of 3° , the mid and lower curves correspond 
to the parent population of quasars with the viewing an- 
gles of 40° and 60°. From the figure, firstly, we find that 
the observed EC component is weakly affected by veloc- 
ity structures compared to synchrotron or SSC ones. The 
reason is that for continual jet the luminosity Lec has the 
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Fig. 10. Under fi apP = 3.6c, the approximately fitting for 
NED data of 0716+714 (Note that we use different param- 
eters for different velocity structures, see the Table 1.) and 
the extrapolated SEDs of parent population. The solid lines 
represent the jet with uniform velocity structures (§ = 0) and 
the dash lines denote the jet with velocity structures (§ = — 1). 
The upper curves report the emission at angles of 2°, the red 
and green curves present parent population of BL Lacs with 
the viewing angles of 40° and 60° . 
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Fig. 11. Under (3 a p P = 3.6c, the approximately fitting for 
NED data of 0716+714 (Note that we use different param- 
eters for different velocity structures, see the Table 1.) and 
the extrapolated SEDs of parent population. The solid lines 
represent the jet with uniform velocity structures (£ = 0) and 
the dot lines denote the jet with velocity structures (£ = 1). 
The upper curves report the emission at angles of 2°, the red 
and green curves are present parent population of BL Lacs 
with the viewing angles of 40° and 60° . 
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Fig. 13. The ratio of luminosity between different viewing 
angles under different £ and bands for 0716+714. 
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Fig. 12. The ratio of luminosity between different viewing 
angles under different £ and bands for Mrk 421. 



,(3°)/L o (40°) M L(3°yL(60 o ) 
L (3>L(40°) ™ A YA 

I I I I I | 

I 



I A J ^ 



Fig. 14. The change of U £ /Ub with z in uniform and deceler- 
ation velocity structures. The dot line is U e = Ub , solid (£ = 0) 
and dash lines (5 = — 1) show Mrk 421, and the dash-dot 
(£ = 0) and dash-dot-dot lines (£ = -1) denote the 0716+714. 

relation of Lec ~ N e 6 3 j EC ~ N e S 3 U' ext ~ T 2 S 3 U ext N e , the 
integration of T 2 5 3 over z gives less discrepancy between 
the models of £ = —0.5 and £ = 0. 

7. Discussions and Conclusions 

In this work we mainly pay attention to the influence 
of the bulk velocity structures on radiative spectra and 
ignore the radiative energy loss of the electrons in the 
studied jet scale. It is reasonable if we suppose ongoing 
particles acceleration in situ can compensate the radiative 
loss. However, radiative losses become important for high 
energy electrons producing the high energy part of syn- 
chrotron and IC emission, when the particles acceleration 
mechanism do not offset the radiative loss. In this situa- 
tion, the high energy part of synchrotron and IC emission 
will be mainly originated at the base of jet, and their spec- 
tra will be steeper with the distance. 

In this work, the synchrotron self-Compton (SSC) mod- 
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Fig. 15. Under (3 a pp = 3.6c, the approximately fitting for 
data of Mrk 421 (from Macomb ct al. (1996) and NED) and 
the extrapolated SEDs of parent population. The solid lines 
represent the emission of jet with uniform velocity structures 
= 0) and the dot lines denote the jet with velocity struc- 
tures (f = 1). The upper curves show the emission at angles 
of 3°, the mid and lower curves present parent population of 
BL Lacs with the viewing angles of 40° and 60°. 



Fig. 17. Under /3 app = 6.8c, the fit curve for data of FSRQ 
3c 279 from the Inoue & Takahara (1996) with the models of 
uniform (£ = 0) or velocity structures (£ = —0.5) (Note that 
wc use different parameters for different velocity structures, 
see the Table 1.) and the extrapolated SED of parent popu- 
lation. The solid lines represent the jet with uniform velocity 
structures (£ = 0) and the dot lines denote the jet with veloc- 
ity structures (£ = —0.5). The upper curves show the emission 
at the angle of 3°, the red and green curves corresponds to 
the parent population of quasar with the viewing angles of 
40° and 60° . 
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Fig. 16. Under (3 a pp = 3.6c, the approximately fitting for 
NED data of 0716+714 and the extrapolated SEDs of par- 
ent population. The solid lines represent the jet with uniform 
velocity structures (f = 0) and the dash lines denote the jet 
with velocity structures (£ = 1). The upper curves show the 
emission at angles of 2°, the mid and lower curves present 
parent population of BL Lacs with the viewing angles of 40° 
and 60°. 



els and External Compton (EC) models of AGN jets 
with continually longitudinal and transverse bulk veloc- 
ity structures are constructed. But our models can not 
discriminate a lateral and longitudinal deceleration sce- 
nario. For the former scenario we might observe the 
limb brightening in VLBI (Giroletti et al. 2004a). For 
the longitudinal deceleration scenario it should exist an 
outward-increasing radio-to- X-ray ratio (Georganopoulos 
& Kazanas 2003b). But they might be interrelated, for ex- 
ample, in a faster spine and slower layer model when the 
IC power is compared with the total bulk kinetic power. 
The Compton rocked effect will cause the spine to recoil 
and decelerate, and form the longitudinal velocity struc- 
tures (Ghisellini et al. 2005). Similarly, if the gradual 
entrainment from external medium cause the jet deceler- 
ation, it would produce velocity gradients across the jet. 
A faster spine and a slower sheath are formed. 

The rapid TeV flares observed for PKS 2155-304 and 
Mrk 501 (Aharonian ct al. 2007; Albert et al. 2007) ap- 
pears to indicate large bulk Lorcntz factors. The system- 
atic differences of the Doppler factors inferred from differ- 
ent waveband observations indicate that the jet holds ve- 
locity structures within sub-parsec to parsec scale. Based 
on our calculation, large bulk Lorentz factor in the base of 
jet for longitudinal velocity structures or in the spine for 
transverse velocity structures, can satisfy the TeV photon 
escape, where the apparent speed is moderate. 

In the unified schemes, we find that beside the view- 
ing angles, the bulk velocity structures play an important 
role in the Doppler beaming pattern. Firstly, from the ob- 
served data to deduce the intrinsic quantities, the velocity 
structures and the viewing angles should be considered at 
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Fig. 18. For different bulk Lorentz factor of 30, 20, 10 and 
5, the Doppler factors vary with the view angles. 

one time. In fact, if the jets have velocity structures, 
their intrinsic flux will be difficultly deduced under dif- 
ferent viewing angles. Secondly, the velocity structures 
also bring some statistical discrepancy between observa- 
tion and theory in simple one-zone model (Chiabcrge et al. 
2000; Henri & Sauge 2006). In one-zone model, if 6 < l/T, 
one has 1 < 8 < 2T except 5 l/T. This means that for a 
few Doppler amplified sources, one expects a large num- 
ber of unbeamed counterparts. In this work, we find that 
under larger viewing angles the velocity structures lim- 
ited by apparent speed have larger effect on the observed 
flux. The statistic of the flux limited sources will provide 
a framework to explore the velocity structures of AGN 
jets. 

In the Fig. 18 also given by Urry & Padovani (1995), 
we present the relations between Doppler factor, Lorentz 
factor and viewing angle. The curve lines for different 
Lorentz factors cross each other within the viewing angles 
of 1° to 10°. They show that large bulk Lorentz factor 
docs not always produce large Doppler factor when the 
viewing angle increases. Furthermore, the Doppler factor 
with large Lorentz factor decreases quickly under the in- 
crease of viewing angles. These properties cause complex 
effects on emissive spectra. 

The observed apparent speed might be the results of 
flux convolution over the resolution region, which can 
solve the inconsistence of bulk motions obtained by differ- 
ent wavebands. We find that the influences of the velocity 
structures upon the observed spectra heavily depend on 
the viewing angles. Under the jet with velocity structure 
constrained by jet power and apparent speed, we show 
that the SEDs have large difference for the jets with and 
without velocity structure in large viewing angle. For the 
jet with transverse and longitudinal deceleration velocity 
structure, its flux is much larger than one with uniform 
structure. Especially, in longitudinal bulk velocity struc- 
tures, as the elementary volume expansion (£ > 0) or com- 
pression (£ < 0) greatly change the electron energy distri- 



bution, we use smaller values of Nq, ^max and B than the 
uniform ones to fit the SEDs; on the contrary, we must 
use larger values for the accelerating models to get fitting. 
We find that for the deceleration velocity structures the 
U e is larger than Ub at the inner of studied jet scale, along 
the jet, Ub is became dominant at the outer. The EC ob- 
served spectra are weakly affected by velocity structures 
compared to synchrotron and SSC ones. 

Nowadays, the VLBA can not resolve the velocity 
structures with sub-parsec to parsec scale, the AGN's 
unified schemes exit a large uncertainty. The TcV 
emission is easily produced in the jet with large velocity 
structures and indicates that TeV blazars should have the 
jet bulk velocity structures within sub-parsec to parsec 
scale. 
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Table 1. Input parameters of the models for the Mrk 421, 
0716+714 and 3c 279. The unit of the zq, Zm ax , and R are 
cm; B is Gauss 
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